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Abstract 

We study the effect of edge diffraction on the semiclassical analysis of two 
dimensional quantum systems by deriving a trace formula which incorporates 
paths hitting any number of vertices embedded in an arbitrary potential. This 
formula is used to study the cardioid billiard, which has a single vertex. The 
formula works well for most of the short orbits we analyzed but fails for a few 
diffractive orbits due to a breakdown in the formalism for certain geometries. 
We extend the symbolic dynamics to account for diffractive orbits and use it to 
show that in the presence of parity symmetry the trace formula decomposes in 
an elegant manner such that for the cardioid billiard the diffractive orbits have 
no effect on the odd spectrum. Including diffractive orbits helps resolve peaks 
in the density of even states but does not appear to affect their positions. An 
analysis of the level statistics shows no significant difference between spectra 
with and without diffraction. 

PACS numbers: 03.20.+i, 03.65.Sq, 0545 
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I. INTRODUCTION 



Periodic orbit theory provides a method of relating local, canonically invariant infor- 
mation about classical periodic orbits to global quantum information such as the density 
of states. However, this theory must be extended if the classical mechanics is not defined 
due to discontinuities. There is one class of discontinuity which is relatively mild in that 
all trajectories are well defined but just their behavior changes abruptly at some points in 
phase space. Examples of this include grazing angles in billiards and the circle-straight 
joint in the Bunimovich stadium @J^. A more severe discontinuity is one in which some 
trajectories are undefined. Examples of this include the vertex of a wedge [^|-^] and three 
body collisions [0 since in neither case can we continue a trajectory through the discon- 
tinuity. Other examples of discontinuities include scattering singularities [|T^], flux tubes 
|l^ , p!3[| and small scattering disks U3|,|l4|. In each case, periodic orbit theory can be ex- 



tended by incorporation of so-called diffractive effects - in the case of vertices this is called 
edge diffraction. To incorporate the effect of a discontinuity, one compares to the solution 



of the local scattering problem. For the wedge this was solved by Sommerfeld |]15[ and the 
solution discussed in Refs. ||T6| , p!7[ . 

The structure of the paper is as follows. In Section II we derive a trace formula for 
diffractive orbits analogous to the Gutzwiller trace formula for ordinary orbits. The ampli- 
tude of each diffractive orbit is affected by the curvatures it experiences on the geometric 
part of its path as well as by the diffraction at the vertex. As an example, we employ the 



theory in a numerical study of the cardioid billiard |]T8[ which is ergodic [0, and in Section 
III we discuss various properties of this billiard. The comparison of the theory and numerics 
takes place in Section IV. A brief analysis of the spectral statistics of the cardiod billiard is 
presented in Section V, while Section VI contains the conclusion. 



II. TRACE FORMULA FOR DIFFRACTIVE ORBITS 

The semiclassical formula for the trace g{E) of the Green function G{E) of a chaotic 
Hamiltonian is |I| 

g{E) ^ TrGiE) = r^^\ ^ eM^{S,/h - a,7r/2)}, (1) 



where the sum is over all periodic orbits 7. The factors Ty, S*^, A^ and are the canonically 
invariant periods, actions, stabilities and Maslov indices of the orbit 7. (The canonical 
invariance of is proved in Ref. |]20|.) The — /+ factor refers to whether the orbit is 
direct hyperbolic or inverse hyperbolic. We extract the density of states p{E) from the trace 
through the identity p{E) = —lm[g{E)]/7r. Periodic orbits are singled out because they have 
a stationary phase with respect to small deviations. The requirement of stationary phase 



can also select other phase space structures and we will henceforth refer to the orbits 



which enter Eq. (|I]) as "geometric orbits" to distinguish from other possibilities. 

Another process which can lead to a stationary phase is a trajectory which hits a vertex. 
While not a classical trajectory, it is still a path in the sense of path integrals with an 
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amplitude that can be found by comparison with the scattering solution of a wedge; we call 
such a path "diffractive" . The asymptotic (in h) contribution to the Green function of the 
Scrodinger equation arising from the path from x' to x via the vertex at ^ is P,p!6|JT7[ 



G,{x, x', E) ^ ^d{9, 9')Gf{x, E)Gf{^, x', E), (2) 

where Gj is the Green function in the absence of the wedge and d{6, 9') is called a diffraction 
constant. It equals 

MO o'\ = _o^M:!L/A[ \ + \ \ 

^ ' ' u \ cos(7r/z/) - cos((^ - e')/p) cos(7r/z/) + cos((0 + 6')/^) J ' ^ ' 

where 9 and 9' are measured with respect to the wedge normal (unlike in Refs. |^^) and 
z/ = a/vr, a being the opening angle of the wedge as shown in Fig. |l|a. The +/— sign refers 
to Neumann/Dirichlet boundary conditions. We will mainly be interested in a = 27r for 
which 

d(9,9')==ec(^^)T=ec(i±^). (4) 

Note that d{—9, 9') = ^d{9, 9') and that (i(0, 0) equals 2 for Dirichlet boundary conditions 
and for Neumann. The factor h'^ /2m appears in Eq. (0) because we are using energy de- 
pendent Green functions; it is absent if we use the Green functions of the Helmholtz equation 
in two dimensions. The Green functions of Eq. (0) all have units l/[Energy\[LengthY'. 

For free motion in two dimensions, Gf{x2,xi,E) = —ih^HQ^\k\x2 — xi|)/8m, where 
iJo^"* is the outgoing Hankel function and k = y/2mE/fi. This form of Gf is assumed in the 
derivation of the diffraction constant. However Eq. (|^) is more general. In the presence of a 
potential with a vertex, we do the local scattering problem by assuming the potential does 
not change much in a typical wavelength. We would then call the directions of the trajectory 
when it enters and leaves the vertex 9' and 9 and use these in determining the diffraction 
constant d{9, 9'). Away from the vertex, we connect the outgoing free space Green functions 
to the relevant semiclassical ones for that potential. For example, Eq. (H) is valid for the 
problem of motion bounded within a wedge in the presence of gravity ||2^ for which G/ is 



more complicated. For that reason, we do not assume billiard conditions in the subsequent 
discussion, although we do restrict the discussion to two spatial dimensions. The content 
of Eqs. (ID and (|^) is that an orbit entering the vertex can be continued out at any angle 
with a quantum amplitude given by d{9,9'). The contribution of such a diffractive orbit is 
of order \/h relative to that of a geometric orbit. We now analyze this semiclassically to 
derive a trace formula in analogy to Eq. (|l|). 

If there are n diffractions as sketched in Fig. Db, Eq. (H) generalises to 



(5) 



where we have suppressed the energy dependence in G j and the 9 dependence in the diffrac- 
tion constants. To obtain the trace, we first identify the points x and x' and then invoke 
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stationary phase to require that the momenta also match smoothly - as for geometric orbits. 
Orbits which satisfy these constraints we call diffractive periodic orbits. They can be found, 
in principle, by firing out trajectories at all angles from all the vertices and determining 
which ones return to a vertex. We impose no constraint on the momenta at the vertices and 
so allow any incoming or outgoing angles. 

To proceed, we define local coordinates along the various classical paths. At each point x, 
we take z and y to be the local coordinates parallel and transverse to the path, respectively. 
At each vertex we define local coordinates Q and rji on the the incoming path and Q and 
ri[ on the outgoing path. The local transverse momentum at a; is p and the local transverse 
momenta at the vertex are vTj and tt^'. The constant energy approximation to the classical 
Green function from any xi to any X2 is fl]] 

Gf{x2, xi, E) ^ -4 X— D{x2, xi) ey.Y>{i{S/h - /i7r/2)}. (6) 
III \/2mh 

The action S is evaluated along the path and /i counts the caustics along the orbit. The 
factor D equals pl 



D{x2,Xi) 



\Z2ZI 
1 



1/2 



dy2dyi 



\z2Z1 



11/2 




1/2 



'1/2 

(7) 



where the subscript on the bracket of the second equation indicates that we take the deriva- 
tive of the final position with respect to the initial momentum while holding the initial 
position fixed. This then defines a fan of initial conditions radiating from the source point 
Xi. This is the contribution of a single classical trajectory - if there is more than one, we 
must sum over them. 

Evaluation of the trace involves integrating along the periodic orbit and transverse to it. 
The integral along the orbit can be done one arc at a time and below we consider just the 
arc between ^1 and ^2- Eqs. (H) and imply 

G;(6,x)G^(x,ei)^ f4^==fl I^(6,a;)/^(x,ei)exp{^[^2i//^-(/i2 + /i2)vr/2]}, (8) 

where 5*21 = 82 + 81 is the action of the path between ^1 and ^2 via x. At each point x along 
the orbit, we calculate the transverse y integral. The only y dependence, to leading order in 
is in the action which we approximate as 



« + (9) 



The partial derivative is taken while holding the initial and final coordinates 77^ and ri2 fixed 
at zero. 821 is the action evaluated at ?/ = and is independent of the position z along the 
orbit. The stationary phase integral yields 
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/ dyexp {iS2i/n) = — |i exp (i{S^Jh - an/2)) , (10) 

where a is zero if the second derivative in Eq. (j^) is positive and is unity if the second 
derivative is negative. 

We now seek to manipulate the various partial derivatives which come from the two 
amplitude factors D and from the stationary phase integral (|10D- These all come with a 
power of —1/2 and with absolute value signs. For purposes of manipulation, we neglect 
those for the moment so the combination we need to analyse is 



dy\ ( dy\ d\Si + S2 





'111 



We have used the fact that the derivative of the action with respect to the final position 
equals the final moment while the derivative with respect to the initial position equals the 
negative of the initial momentum. Also, we have inserted all of the relevant subscripts 
to indicate what is being kept fixed in each derivative. We now show that this factor is 
independent of position along the orbit. We can combine partial derivatives in Eq. (|TT|) to 
obtain 



f 9p\ ( dy\ ( dp\ I d{y,p) 



dTT2 I \ On', I , \ dn', I , \ dn2 I \ d(n2, it',) I , 



(12) 



In the right half of Eq. ([T2|) we have borrowed the Jacobian notation of Ref. ||2^. (We have 
made use of the trivial freedom to specify that in the first derivative of the left hand side 
we are also holding rj^ fixed, with similar specifications in all four terms.) To determine 
the value of this Jacobian corresponding to some different point z' along the trajectory, we 
should multiply Eq. (^) by the Jacobian relating the transverse variables variables (|/, p) at 
position z to the transverse variables {y',p') position z'. However, because these variables 
are transverse to a trajectory, their Jacobian is identically unity; the two sets of variables 
are canonically related owing to symplectic nature of the Hamiltonian flow. It follows that 
the combination of factors appearing in Eq.([TT|) is independent of position z'. In particular, 
it is particularly convenient to calculate it very close to one of the vertices. If z is such that 
the point is close to ^2, "we have 



It follows that 

/ d{y,p) 



,5(7r2,7rl) 




(14) 



where the right hand side is evaluated at the point ^2- Heceforth, we change notation slightly 
and call this term |p- to stress that it is evaluated at the second vertex. This factor is simply 
the spread in position at .^2 of a fan of trajectories radiating from ^1. 
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We also want that the phase index /ii + /i2 + o" be independent of position along the orbit. 
It is not true that the indices are separately invariant; it is simple to imagine that as we 
change position z along the orbit, we will gain or lose caustics in going from the two vertices 
to the intermediate position. However, these changes will be exactly mirrored by changes 
in the index a such that the sum is invariant. For a demonstration of this, we refer to Ref. 
| 25| where the authors evaluate an integral similar to Eq. (|^). They interpret Eq. @ as a 
propagator along the orbit with z playing the role of time. Using the semigroup property of 
time-dependent propagators, they conclude that the phase index is a constant. Since this is 
constant, we are free to use [i2\ which is the number of caustics of a fan of trajectories going 
from vertex 1 to vertex 2. 

That completes the y integral. For the z integral, we remark that the only z dependence 
is in the velocity which appears in the amplitudes D of Eq. (^). Since there are two 

of them, the integral to be performed is simply 



^ =T2i 



(15) 



which is just the time it takes to get from to .^2- Putting together all the remaining 
factors, we conclude that 



y dzdyGf{^2,x)Gf{x,^i 



T, 



21 



(z^)2v/M|C2Cl|'/' 

%GK6,ei), 



dr]2 



dTi[ 



-1/2 



exp{z(S'2i - /i2ivr/2)} 



(16) 



where we have used Eqs. (^) and (|^) in the second line and dropped the superscript on 
the action. 

The appealing fact that the trace integral on the arc between and ^2 is proportional 
to the Green function between these points simplifies the analysis tremendously. Recall that 
we must still multiply all of the other free Green functions from Eq. (^ so there is a factor 
which is simply the product of all the Green functions from vertex to vertex. In doing the 
integrals along the arc between and ^j+i, we get exactly the same product but multiplied 
by Tj+i^j so that the integral of Eq. (1^) is 



2m J 

nj/2 



ih \ Sum? 




drj, 



i+i 



dn' 



-1/2' 



exp {i{S-y/h — ay'K/2 — Sn^ir/i)} , (17) 



where S-y, and are the sums of Si+i^i, Tj+i^j and /ij+i,i along the orbit. The velocity 
(i is given by energy conservation and is a constant at each vertex and the index i is cyclic 
so vertex + 1 is identified with vertex 1. This diffractive trace formula is the main result 
of this section. The formula was given in Ref. by comparison with creeping diffraction, 
where Watson contour integration can be used to show that the trace has the same 
structure. The diffractive trace formula is similar to the Gutzwiller trace formula Eq. (|l]) 
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but is suppressed by a relative factor of h^^"^ Eq. (^) only shows the contribution of 
a single diffractive orbit; in practice we must sum over all such orbits and so introduce the 
subscript 7 to distinguish them. If the orbit is a repeat of a shorter primitive orbit, the 
factor of T in Eq. ( |17| ) is the period of the primitive orbit. 

We now specialise to the potential-free case so that h = m = 1, S/h = kL, T = L/k, 
E = k'^/2, and = k. We further invoke the infinitesimal relation Stt'^ = k6(j)'^ where 6(j)[ is 
the angular deviation from the periodic orbit on leaving vertex C,i so that 



dfji+i 


1 




d< 


~ k 





(18) 



The contribution to the density of states in k is given by p^yik) = —klm[g^{E)]/n so 
to be contrasted with 

r 1 

. cos (kLg - agn/2) (20) 



TT 



A, T 1/jAs 



for geometric orbits. In Eq. (|T9D the factor Fi has a simple interpretation; if we launch a 
narrow cone of trajectories from vertex centered on the periodic orbit, Fi gives the width 
of the cone when it arrives at C,i+i [0. This interpretation in terms of cones is in contrast 
with that of cylinders for geometric orbits. Eq. (|19D was also obtained in Ref. for the 
special case of straight walls everywhere so that Fi = Li, the distance between the vertices. 
We also mention that this analysis applies equally well if two of the diffraction points are 
at the same vertex. In particular, if there is only one vertex, then the diffractive periodic 
orbits are those which leave the vertex and return to it following a classical path. Eqs. (|1^ 
and (|19|) are true regardless of whether the classical motion is chaotic or not, although they 
are restricted to isolated diffractive orbits. 

There is a zeta function p6| corresponding to Eq. (0) in analogy to that which exists for 
geometric orbits [^. We start by defining the diffractive weight for each diffractive orbit 



n-y/2 




drj. 



'i+i 



-1/2' 



exp {i{Sy/h 



/2-3n^7r/4)} 



(21) 



where we have now introduced the subscript 7 to distinguish diffractive orbits. The sum 
over diffractive orbits is then 



9d{E) 



T 



T 



ih I — u 



(22) 



where we have organised the sum into the primitive orbits and their repeats. To leading 
order in h, 



so that 

9<iiE) = ^\og(lli^-t,)y (24) 

The quantity (^^-^{E) = Yi-yi^ ~ ^7) is the diffractive zeta function to the power —1. When 
multiphed by the corresponding geometric zeta function to the power —1 and appropri- 
ately regulated |]28| . |29[| , the product equals the spectral determinant Yln{E — En) so that its 
zeros give the quantum energy levels. Due to the regularisation, the zeros of the separate 
terms in (^^{E) are not true zeros of the product. A formula analogous to Eq. ( p4D also 
holds in the case of billiards . 

The function (^^(E) involves only a single product over periodic orbits. In contrast, the 
zeta function for geometric orbits has an additional product over an integer index which 
can be thought of as labeling local eigenstates transverse to each orbit [0. Near the orbit, 
these local eigenstates typically have the form ipniu) ~ 2/" with n > 0. We conclude that 
diffractive orbits have only the n = local eigenstate. Higher states do not exist because 
they would have a node on the periodic orbit and would not be affected by the diffraction. 
This was also noted in the scattering geometries discussed in Refs. |^,|^ where it caused 
there to be no lower families of quantum resonances. This difference is intimately related to 
the fact that the Green functions in Eq. (^) are multiplicative |^ in contrast to the behavior 
of Green functions for geometric orbits 



III. THE CARDIOID BILLIARD 

In this section, we discuss various aspects of the cardioid billiard which are relevant to 
us. We briefly review its classical properties and construct the symbolic dynamics for the 
geometric orbits and diffractive orbits. The Maslov indices are given by a simple rule in 
terms of the symbolical dynamics. We then discuss the role of symmetry in the quantum 
problem and how the geometric and diffractive trace formulas conspire to give the even and 
odd spectra. We conclude with a discussion of the Weyl formula. 



A. Classical Mechanics 

We study the cardioid billiard whose boundary is defined by the following mapping of 
the unit circle in the complex plane 

z(^) = e^^ + V^ ^G[-7r,7r), (25) 

and is plotted in Fig. ^. The angle 6 is defined such that it changes discontinuously from vr 
to — TT at the cusp. If the factor of 1/2 is replaced by a parameter b this represents a family 
of billiards introduced in Ref. |18[ and subsequently studied exhaustively p9|j3T| -p5|. For 
b < 1/2, Eq. (p5D is a conformal mapping but is not strictly conformal for 6=1/2 since the 
derivative of z with respect to 6 vanishes at ^ = vr. In practice, this does not matter and 



the algorithm introduced in [33| to find the quantum eigenvalues still applies and was used 
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by us. For 6 > 1/2 the curve crosses itself near 6 = n and the bilhard is not well defined. 
We note that the quantum behaviour of the cardioid was recently studied in Ref. [^]. The 
dynamics in the billiard consists of free motion within the domain followed by equal angle 
(specular) refiections at the boundary. Trajectories which strike the vertex are not defined 
but these are of measure zero. Motion in the cardioid has been proven by Markarian to be 



ergodic |jT9|. It is similar to the Bunimovich stadium in that it is defocusing. Defocusing 
means that each point on the billiard has positive curvature so that parallel rays striking the 
boundary are initially focused. However, the billiard geometry is such that the trajectories 
typically diverge even more after the focal point resulting in a net defocusing. It is this 
mechanism which leads to the average divergence of trajectories and to chaos. In contrast, 
a dispersing system such as the Sinai billiard has negative curvature so that initially 
parallel rays striking the boundary are immediately dispersed. 
The curvature of the cardioid is 

K{e) = lsec(^^y (26) 

which is positive for all 6, as in a circle. For b somewhat less than 1/2, the region near 9 = ir 
is dispersing rather than defocusing. Billiards with mixed focusing properties like this are 
difficult to analyze mathematically [|T^,§^ and it is for this reason that Markarian's proof 
works only for the cardioid. In cartesian coordinates near the cusp the billiard boundary 
satisfies the equation 

y^±^(-(2x + l)f\ (27) 

Consequently there is a cusp at x = —1/2 which locally looks like a half plane extending 
to the left. This is an example of a vertex singularity so that the analysis of the previous 
section applies. 

We found periodic orbits numerically by using the principle of least action. For an 
arbitrary periodic orbit the number of intersections with the boundary was specified and 
the intersection positions were varied until a local minima of the total orbit length was 
found. Diffractive orbits were found the same way but with the constraint that one of the 
intersections was at the cusp. Various geometric orbits are shown in Fig. ^ The label of 
each orbit includes the number of intersections and also a letter index to further distinguish 
them. We describe a better naming system below. The asterix designates self-dual orbits as 
defined below. In Fig. ^ we show various diffractive orbits. The naming scheme is similar 
to before, the number gives the number of bounces - not counting the cusp. The primes 
indicate diffractions, ie. the number of encounters at the cusp. Orbits 2a', 3a" and *4a" 
have arrows drawn to indicate the scattering directions at the cusp. The last two have two 
diffractions and 3a" is seen to be a composition of la' and 2a' 

The three orbits *6b, *8b, and *10b are geometric and reflect specularly near the cusp, 
contrary to appearances. Also, 4a misses the cusp and is geometric in contrast to 4a' which 
hits the cusp and is diffractive. In Fig. ^ we show examples of pruned orbits ||D| of both 



the geometric and the diffractive kind. These have an index p in their label to indicate that 
they are pruned. It can be seen that the orbits are related pairwise; orbit lOp' appears to 
be composed from 5p and 5a' and similarly for 12p' and 14p'. This is a feature we discuss 
below. 
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B. Symbolic Dynamics 



Symbolic dynamics is the partitioning and labeling of topologically distinct regions 
of phase space. Because of the reflection symmetry of the problem, we can discuss the 
dynamics either in the full domain or in just half of it. The half domain, also called the 
fundamental domain, has dynamics which are the same as in the full domain but with a 
reflection at the symmetry axis. We will show that the full and fundamental domains have 
distinct but closely related symbolic dynamics. We begin with a discussion of the symbolic 
dynamics of the geometric orbits. 

For an arbitrary, time-reversal invariant system with a reflection symmetry, all orbits 
belong to one of five classes. In principle there can be boundary orbits which lie directly on 
the symmetry axis; in this example there happen to be no such geometrical orbits. The other 
four possibilities are (i) symmetric and self retracing, (ii) symmetric but not self-retracing, 
(iii) self-retracing but not symmetric, and (iv) neither self-retracing nor symmetric. These 
occur with multiplicities 1, 2, 2 and 4 respectively and examples are *4b, 3a, 6c, and 7b. 
Although (iv) is the most generic possibility, we did not find many examples of it among 
the shortest orbits. It is typical that the shortest periodic orbits are special |]^J^ and 
that generic ones begin to appear only for longer lengths. Orbits of classes (ii)-(iv) behave 
identically in the half domain as in the full domain. Orbits of class (i), the so-called self-dual 
orbits, must be treated more carefully because in the fundamental domain they are periodic 
in T^/2 as well as in T^. 

We start by discussing the symbolic dynamics of the full domain. Recalling that each 
point on the boundary of the billiard is labeled by an angle 9 in Eq. (|25|) , we assign a 
trajectory the symbol "+" every time it has a reflection which increases 6 (counterclockwise) 
and the symbol "— " every time it has a reflection which decreases 6 (clockwise). At the 
cusp the angle 6 changes discontinuously by 27r, so it defines what we mean by an increase 
or decrease of angle. Since we allow no geometric orbit to hit the cusp, it is not a problem 
that the sign of 9 is not defined there. This is a general property of dynamical systems - the 



symbolic dynamics is often conveniently described with reference to a discontinuity |2^ . 
As an example, consider orbit 3a. When going in the counterclockwise sense, it has the 

symbol sequence H h. Its time reversed partner, which is distinct in the full domain, has 

the sequence — I — . We are free to start counting symbols anywhere on the orbit so that 
any cyclic permutation of a symbol sequence describes the same orbit and is not distinct. If 
W is the symbol sequence of an orbit which is symmetric under reflections (such as 3a) then 

W = W*, where W* is obtained by reversing the order of the symbols. For example, H h 

reads the same left to right as right to left. Self- retracing orbits of length 2n (they must 
be even) have symbol sequences with the structure W = AA where A is some sequence of 
length n and A is obtained from A by reversing the order and every sign. For example, the 

symbol sequence for the self-retracing orbit 6c is 1 h-l- for which A = h and 

A = — h +. A self-dual orbit is both symmetric and self-retracing and therefore its symbol 

has both properties. An example is *4b whose symbol sequence is + H (this satisfies 

the property of being symmetric if one makes use of the freedom to cyclically permute the 
symbols.) 

The symbolic dynamics in the fundamental domain are defined by looking at the symbol 
sequence of an orbit in the full domain and assigning a if two adjacent symbols are the 



10 



same and a 1 if they are not. For example, orbit 3a, which is labeled H h in the full domain 

is Oil in the fundamental domain. The time reversed orbit in the full domain, — I — , is also 
Oil in the fundamental domain. This is consistent since in the full domain they are distinct 
and should have separate symbols while in the fundamental domain they are not distinct 
and should have the same symbol. 

Note that a self-dual orbit in the full domain has a symbol sequence in the fundamental 
domain which repeats itself, for example *4b has the sequence 0101 = (01)^. Therefore 
a self-dual orbit, when mapped onto the fundamental domain, is the double repeat of a 
shorter orbit. Every odd multiple of this shorter orbit will be present in the fundamental 
domain but not in the full domain and this is apparent in the symbolic sequences. Any 
sequence in which 1 appears an odd number of times can not be periodic in the full domain 
and so corresponds to a self dual orbit. Any orbit which is symmetric in the full domain is 
self-retracing in the half domain. 

We show in Table I, the symbols of all orbits up to length 4 as measured in the funda- 
mental domain. There is no fundamental orbit 0, this is reminiscent of the co-linear helium 
problem |]10[, among others. The pruned family 5p,6p,7p... means that there are no orbits 
with the symbol sequence 0"11 with n > 2. On the other hand, there is an accumulation 
of whispering-gallery-like orbits labeled 5a,6a,7a... whose symbol sequences have the form 
0"101 and whose lengths accumulate to L = 12 as n — > oo. There are also orbits of the form 
01"; for n even they are 3a,5b,7c... while for n odd they are self-dual and are 4b, 8b,... We 
believe that both series exist for any n. 

The symbolic dynamics of diffractive orbits is clearest if one think of the cusp as not 
being part of the boundary but rather being a means of getting from one point on the 
boundary to another. We introduce a symbol d which represents a path between two points 
on the boundary which goes via the cusp and keep + and — as defined above. Therefore, 
orbit 2a' has two boundary intersections and has a symbol sequence d+ or rf— depending 
on the sense of the rotation. Symmetric diffractive orbits labeled dW satisfy W = W* while 
self retracing orbits have the structure W = AA, as for geometric orbits. 

The rule for the symbolic dynamics in the fundamental domain is again found by looking 
at the word in the full domain. A symbol d maps to a d, while a + or — maps to a or 
1 depending on whether the next non-d symbol is the same or different. For example, in 
the half domain, 2a' has the symbol dO. In the fundamental domain the two rotation senses 
of 2a' are not distinct and it is consistent that there is only one symbol sequence. The 
geometric orbits discussed above can be subsumed into a larger ternary alphabet in which 
they are the subclass with no d in their symbol sequence. 

As before, any symbol sequence with an odd number of I's must correspond to half of 
a self-dual orbit. However, self-dual orbits with one diffraction have the special property of 
being geometrically identical to a non-self-dual diffractive orbit and we refer to the pair as 
complements. An example of this is the self-dual orbit *4a" which is a perfect overlap of orbit 
2a', the only difference being that the first backscatters at the cusp while the second forward 
scatters. Consequently, *4a" has twice the length of 2a' and suffers two diffractions rather 
than one. This is general, the self-dual orbit always leaves the cusp at an angle which is the 
negative of its complement and then follows a trajectory which is simply the reflection of its 
complement. If an orbit in the full domain has the symbol sequence dW then its self dual 
complement has the symbol sequence dWdW where W is defined such that WW satisfies 
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the self-dual property. For example, *4a" has the symbol sequence d + d— and H — is clearly 
self-dual. In the fundamental domain the symbol sequence of a self-dual orbit is found from 
its complement by switching the symbol immediately before the d. For example, the symbol 
sequence of 3b' is dll while that of its self dual complement (not shown) is (ilOrflO (recall 
that due to the cyclic symmetry, the last character in these sequences is "before" the d). 
In Table II we show the labels of all singly diffractive orbits up to length 4. Note that the 
symbols dlOl and ciOll are simply time reversed copies of each other and contribute equally 
to the trace. 

To find the symbol sequence of any multiply diffractive orbit, we use the same rule. For 
example, the doubly diffractive orbit 3a" starts at the cusp, travels along la' back to the 
cusp, diffracts onto the 2a' orbit, travels along 2a' back to the cusp and finally diffracts onto 
the la' orbit in the original direction. Its symbol sequence is d^d. Because there is only one 
vertex, the only multiply diffractive orbits are compositions of singly diffractive ones. More 
possibilities would exist if there were more than one vertex. 

If all possible symbol sequences were realised as orbits, the number of singly diffractive 
orbits in the fundamental domain up to length n would grow as 2". For geometric orbits with 
a complete binary alphabet, the number grows as 2^ /n. The factor of n in the denominator 
is because cyclically permuted symbols correspond to the same orbit and should be counted 
only once. It appears that there are more singly diffractive orbits of long length and it is not 
clear if they ultimately dominate the spectrum. A related issue is whether the pruning of 
geometric and diffractive orbits is such that the exponential proliferation of the two classes 
of orbits is given by the same exponent. This is quite likely since very long orbits cover 
the phase space uniformly and are therefore susceptible to the same pruning mechanisms as 
shown in Fig. ^. These questions will be studied in greater detail in a later publication |43 
and here we have just a brief discussion. 

The pruning of geometric and diffractive orbits appear to be strongly correlated. For 
example, the pruned diffractive orbit lOp' looks as if it is composed of the pruned orbit 5p 
and the diffractive orbit 5a'. This is confirmed by looking at the symbol sequences of these 

three orbits. Orbit lOp' has the word (i + + + H 1- + ++ which is equal to the composition 

of those for 5a' and 5p which are d + + + + and — h + + +, respectively. A similar result 
holds for orbits 12p' and 14p'. Usually, the existence of pruning implies problems in the 
cycle expansion of the zeta function due to the breakdown of shadowing |2^ so the fact that 
orbits and their shadows have disappeared together might prove very useful. An example of 

non pruned shadowing are 56 2a + 3a (since H |-H — = H — || + H — ) and 4c' -v^ 2a' + 2a 

(since d H h = d +\ \ — h). 

The symbolic dynamics was useful to us in guessing the topology of a few missing orbits. 
However, we did not make extensive use of it as we were content to know the shortest orbits 
and these can be found by trial and error. It is useful, however, in the following discussion 
of Maslov indices and symmetry reduction. In addition, a project which required knowing 
many orbits, such as attempting to find semiclassical approximations to many quantum 
eigenvalues, would need to make extensive use of the symbolic dynamics. 
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C. Maslov Indices 



In this subsection we discuss the Maslov indices which appear in the trace formulas 
(H) and (|l3), beginning with the geometric indices. Starting at an arbitrary point x on 
the orbit, the index equals the number of caustics fi plus an index u which arises on 
doing the stationary phase integral in the determination of the trace. Although fi and u 
depend separately on the point x along the orbit, their sum does not. In fact, is a 
canonical invariant which equals twice the number of times that the stable and unstable 
manifolds wind in completing one circuit of the periodic orbit. It follows that we can use 
any point x on the orbit to calculate a^. 

We numerically propagated the 2x2 monodromy matrix M for each geometric orbit 
and counted the number of caustics n by the number of times that one of the off diagonal 
elements M12 changed sign. At the end, we also found the value of (TrM — 2)/Mi2; if it is 
positive u = and if it is negative u = 1. Doing so, we found the following simple topological 
rule; equals the number of reflections or, equivalently, the length of the symbol sequence 
of that orbit. This rule probably arises from the fact that the cardioid is a purely defocusing 
billiard so there is, on average, one focus per reflection. In the Bunimovich stadium, one 
finds a similar rule that is incremented by one for each reflection off the defocusing end 
caps . It is common that such a simple rule exists and it is usually related in some simple 
way to the symbolic dynamics. 

For the diffractive orbits is the number of caustics between successive diffractions. 
We found that always equals the number of geometric reflections and hence the length 
of the symbol sequence, as before. This appealing result implies a unity between the two 
classes of orbits. 

Every geometric reflection also induces a sign change due to Dirichlet boundary con- 
ditions. We account for this by incrementing by 2 at every reflection, so that in total 

= 3m^ where m-y is the number of reflections experienced by orbit 7. For Neumann 
boundary conditions this is not necessary and = m^. Consistent results, specialised to 
the half cardioid, were obtained in Ref. p6|. 



D. Symmetry Decomposition 

The billiard has a reflection symmetry C2 and consequently all quantum states can be 
classified as even or odd. The trace decomposes as 

g{E) = g^{E)+g,{E). (28) 

We can separately find g±{E) by studying the dynamics in the fundamental domain and from 
them the densities p±{E). The behavior of the geometric trace formula (|lj) under symmetry 
decomposition is a well studied problem pT|^^^ and here we review the results which 
are relevant to us. Thereafter, we discuss the decomposition of the diffractive trace formula, 
which is slightly different. 

The non-self-dual orbits have multiplicities of either 2 or 4 in the full domain and their 
amplitudes are divided equally between the two parity classes. Self-dual orbits require more 
care. Half of such an orbit, being periodic in the fundamental domain, contributes to the the 
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separate traces as follows. Its period, stability and Maslov index are all half of the full orbit 
and its stability is the square root. In addition the =p factor in the denominator of Eq. (|l]) 
is replaced by ±. Finally, there is a group theoretical weight of ± corresponding to the 
even/odd parity. This last factor ensures that the contribution of this half orbit identically 
cancels when we evaluate the sum g{E) = g+{E) + g^{E). This is consistent since the half 
orbit is not a periodic orbit of the full domain and should not affect the total density of 
states. The double repeat of a half orbit is the full orbit and its amplitude is divided equally 
between the two parities. 

We next discuss what happens to the diffractive orbits in the presence of the symmetry. 
There is a diffractive boundary orbit; it contributes only to the even spectrum [[7|,P,^ 



unlike a geometrical boundary orbit which contributes to both. The distinction can be 
traced to the fact that the diffractive Green function (|^) is multiplicative in the direct 
Green functions. For the other diffractive orbits, we recall the previous discussion that each 
of them has a self dual complement of twice the length. The only difference between the 
orbits in the fundamental domain is that they have different diffraction constants. If one 
has the diffraction constant d{6,6') the other has the diffraction constant d{—6,6'). Since 
they are in all other respects identical, we can include both of them by defining separate 
diffraction constants for the even and odd parities, namely 

d±{e,e') = d{e,e')±d{-e,0'). (29) 

The different sign for the two cases is the same group theoretical weight mentioned for 
self-dual geometric orbits. 

As mentioned above, for the half plane and for Dirichlet boundary conditions, d{—6, 9') = 
d{6, 6') so we have 

d+{e, 6') = 2d{e, 6') d.{e, e') = o. (30) 

This implies that the odd spectrum is completely insensitive to the existence of diffraction. 
For wedges which are not half-planes, Eq. ( pOf ) is not true but it will still be true that Eq. ( pQ]) 
will cause the two parities to be affected differently. If we study the billiard with Neumann 
boundary conditions we reach the opposite conclusion. First, the boundary orbit has zero 
amplitude since Eq. (H) implies that d{0, 0) = in that case. Also, d{—9, 9') = —d{9, 9') for 
Neumann boundary conditions so that 

d+{9, 9') = d^{9, 9') = 2d{9, 9'). (31) 

Therefore, the odd spectrum would have diffractive peaks and the even spectrum would not. 



E. Weyl Formula 

The Schrodinger equation for a billiard reduces to the Helmholtz equation 

(V^ + A;2) V^(r) = (32) 

with Dirichlet, Neumann or mixed boundary conditions on the boundary of the domain. 
Finding the eigenvalues k"^ leads to the density of states p{k) = J2n ^ik — kn). One commonly 
decomposes this into a smooth part and a fluctuating part 
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p{k) = p{k) + PQ{k). 



(33) 



These terms have distinct classical interpretations. The first term, commonly called the 
Weyl term, is related to the geometry of phase space, such as the area, perimeter, curvature 
and other properties of the billiard boundary. The second term is given by the dynamics as 
encoded in the trace formulae (|I|) and ([T7|). Actually, each term of Eq. (|3^) is an asymptotic 
expansion in powers of [0. To date, the first 16 terms of the expansion of p{k) have 
been calculated but here we content ourselves with the first three. Additionally, the first 
corrections to Pfi{k) have also been determined but we do not consider them here. 

Instead of p{k), one often refers to the spectral staircase function N{k) of which p{k) is 
the derivative. Its expansion for Dirichlet boundary conditions is (see for example Ref. p7| ) 

N(k) ^ —e -—k + C (34) 

^ ^ Arc Att ^ ^ 

A is the area of the billiard, L is the length of the perimeter, and C is related to the curvature 
and to corners by 



c=^i:^+^y«(^M». (35) 

The sum is over angles in the billiard boundary; we have one angle of 2tt. The integral gives 
the total curvature over the boundary of the billiard. Note that Eq. ( P^ also applies to 
billiards with Neumann boundary conditions if we multiply every other term by —1. For 
the full spectrum, we find A = 3Tr/2, L = 8 and C = 3/16. We can also use Eq. (|^) for 
the odd spectrum by taking the billiard domain to be the fundamental domain. We then 
have A = Sn/A, L = 6 and C = 3/16. The difference between the total spectrum and the 
odd spectrum must correspond to the even spectrum, so for it we have A = Sn/A, L = 2 
and C = 0. The symmetry decomposition of the Weyl formula for billiards was discussed 
in much greater generality in Ref. ||50|. These results for the cardioid billiard were also 
obtained in Ref. 

Comparing the exact staircase functions with the approximation (Q) is a useful check 



2 e? 



that there are no missing levels. In addition, the point where there the curves start to 
deviate is a useful criterion for establishing when the numerical eigenvalues are no longer 
reliable. The Weyl formula is also needed to compare the Fourier transforms of the data 
and the trace formulae in Section IV and also to renormalise the spectrum for the statistical 
analysis of Section V. 



IV. NUMERICAL RESULTS 

In this section we present comparison between the exact spectra and the results of pe- 
riodic orbit theory. We first do this by directly comparing the Fourier transforms of the 
exact spectra and the trace formulae in the reciprocal space of orbit lengths L. Overall, 
there is good agreement and we successfully reproduce geometric, diffractive and doubly 
diffractive peaks as well as the interferences among them. However, there is a region of L 
which is not well reproduced for reasons we explain. We also find that for other regions of 
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L, the exact diffractive peaks have magnitudes larger by a few percent than what we expect. 
We explain this by recalling that the cusp is not a perfect vertex. Finally, we use the sum 
over periodic orbits to find the first few quantum eigenvalues and find the odd result that 
including diffractive orbits seems to shift them only slightly. 



A. Fourier Transforms 



Using the algorithm of Ref. we truncate the Hilbert space to contain only the lowest 
6600 states and calculate the lowest 1000 eigenstates of each parity with an accuracy better 
than 0.001 times the average level spacing. It is these 2000 states we use in the analysis. 
For a precise comparison, it is best to work with Pf}(fc) which is obtained from subtracting 
the Weyl term N{k) of Eq. ( P^ ) from the exact density of states. We obtain its Fourier 
transform as 

F{L)= dkw{k)[p{k)-p{k)]e'^^ (36) 



where w{k) is a window function. We chose to use the 3-term Blackman-Harris window ||51 
which gives a good compromise between narrowness of peaks and smallness of side lobes. 
This is defined as 



w{k) 



where (oq, Oi, 02) = (0.42323, —0.49755, 0.07922). This function goes smoothly to zero at /cq 
and ki which we choose as the of the first and last eigenvalues in our spectrum. 

We apply the same Fourier transform to the trace formulae Eqs. (|19D and ( ^OD to obtain 
the semiclassical approximation Fqq[L). We included all relevant period orbits with L < 11; 
the most was for the even spectrum which had 38 orbits, including halves of self-dual orbits 
and multiple repeats of 2a. The results up to L = 10.7 are shown in Fig. ^for the even, odd 
and combined spectra. Other than the region around L = 7.5 which we discuss later, the 
agreement is very good. This was also observed in Ref. P| but here we are also verifying 
the geometrical factors F in Eq. (p!9|). Note that there are no diffractive peaks in the odd 
spectrum as we argued above. The geometric peaks in the even and odd spectra near 
L = 2.5, L = 4.7, ■ ■ ■ are halves of self-dual orbits and are absent in the full spectrum. 

We stress that the relative heights of the diffractive peaks are artifacts of the range of the 
quantum spectrum we choose to consider and should not be used to estimate the relative 
weight of these orbits in determining quantum eigenvalues. The reason is that they are 
suppressed as 1/ \/k and contribute more and more weakly to the energetic states. However, 
their effect in the ground state region can be quite large - a possibility we explore later. 

We now turn our attention to the large discrepancy near L = 7.5. There are two 
distinct reasons why the trace formulas fail there. The first is that orbits 4a, 4a' and 
*10b are close in configuration space and interfere differently than we have assumed up to 
now. The second reason is that the cusp is only approximately a half-plane vertex. To 
understand the first point, we appeal to the calculation which gives the diffraction constant. 



Sommerfeld [15| showed that an incoming plane wave is broken up by a half plane vertex 
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into three components; a plane wave which continues in the original direction, another plane 
wave coming from reflection off one face, and a third component which he identified as the 
diffracted field. The diffracted field is asymptotically an outgoing circular wave of the form 
f {6) exp(ikr) / ^/kr and the diffraction constant d is proportional to f{9). However, in the 
directions close to the two out going plane waves, the diffracted wave takes longer and longer 
to obtain its asymptotic form and exactly in the plane wave directions it never does. It is 
this attempt to connect to the incorrect asymptotic form in these two directions which leads 
to the divergence of Eq. when 6 ± 6' = tt. Ideally, one would like to have a uniform 
approximation to cover all ranges of 9 and 9'. In the case discussed here, we will find that 
the approximation improves as k increases. 

The second reason for the failure of the trace formula near L = 7.5 is that the billiard 
domain departs relatively rapidly from its local half plane geometry at the cusp so that 
the diffraction constants derived assuming a half plane may not be appropriate. In fact, 
this approximation is fine for orbits which come in from the right, such as la', since they 
are never close to the faces of the cusp. However, orbits entering the cusp from the left, 
such as 5a', are very sensitive to this approximation. The curvature of the boundary means 
that such orbits are not as far from the boundary as is assumed in the calculation of the 
diffraction constant so that the Dirichlet boundary conditions cause more suppression than 
is accounted for. 

To understand this another way, note that for any finite wave number k there is not 
infinite spatial resolution so the cusp appears as a finite angled wedge. However, a finite 
angled wedge has different diffraction constants. Therefore, it is reasonable to suppose that 
at small wavenumbers the lack of resolution, inherent in the finiteness of k, manifests itself 
as inaccuracies of peak amplitudes in the Fourier spectra. We explored this effect for orbit 
la' by looking at two windows of k. For this study, it is best to use the combined spectrum 
because then there is no interference from the self-dual orbits 2b and 4b. We are extremely 
sensitive to errors because we are interested in the difference between two small peaks. 
Accordingly, we use just the lowest 285 states which are extremely accurate (approximately 
140 of each parity) which we divided into two windows of approximately equal extent in k 
(2.01-26.32 and 13.09-39.39 respectively with means (k) = 14.17 and 26.24.) The comparison 
between the numerics and the trace formula |l^ is shown in Fig. |^. In both cases, we find 
that the exact spectrum has a peak which is slightly larger than predicted but that the 
discrepancy is larger by about 50% for the first window. 

This strengthens our argument that the large discrepancies around L = 7.5 in Fig. |^ 
arise because the orbits responsible for those peaks are sensitive to the fact that the cusp is 
not a true half plane. We argue that because they approach the cusp from the left they are 
far more sensitive to this than orbits approaching from the right. Furthermore, such orbits 
are probably sensitive to the entire curved geometry in the neighborhood of the cusp. 

The fact that there is structure in the odd spectrum with lengths the same as for diffrac- 
tive orbits (particularly visible near L ^ 7.5, an effect also observed in Ref. 0]) is not 
unexpected. In one of their seminal papers, Balian and Bloch [^] argued that for a billiard 
domain, with a discontinuity in the n'th derivative of the boundary, there exists contributions 
to the density of states of order h"'^'^. For sharp corners, we have n = 1 and this conforms 
to the previous discussion. For discontinuous changes in the curvature we have n = 2. The 
second case applies here for the odd spectrum since the half cardioid has continuous slope 
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but its curvature changes discontinuously from zero to infinity at the vertex. Normally, this 
structure is too small to be visible but we suppose that the geometrical considerations which 
amplify the diffractive peaks in the even spectrum have a similar effect on the higher order 
peaks in the odd spectrum. 

Another possibility for contributions due to a change in curvature is discussed in 



Refs. IQJ^ in the context of the Bunimovich stadium It was found that although 

one can unambiguously continue a classical trajectory which encounters the discontinuity, 
the stationary phase integral to evaluate its contribution to p{E) requires more care. This 
is because the character of the motion is different on each side of the orbit. This applies to 
orbits which reflect specularly at the vertex and so are geometrical - not diffractive. 

Finally, we discuss the small peak near L = 9.8 corresponding to orbit 3a". We use 
Eq. (|1^) with the appropriate diffraction constants. To eliminate interference from the 
nearby self-dual peak 8c we analyze the combined spectrum. Because this is such a small 
peak, it is particularly sensitive to small backgrounds. One source of background is slight 
errors in the determination of the exact quantum energy levels; such errors manifest them- 
selves as weak oscillations in the Fourier transforms. For this reason, we used only the lowest 
1000 levels. In Fig. |], we show both the absolute value of the Fourier transform and its real 
part. The discrepancy, which grows towards the right of the figure, comes from a slight error 
in the peak amplitude of 4b'. 



B. Recovering Quantum Eigenvalues 

The trace formulae Eqs. (|T^ and pOj ) can be used to find quantum eigenvalues. Here 
we do this with and without the diffractive orbits to see what effect they have on the 
determination of the eigenvalues. We simply summed over the orbits used to obtain Fig. |^. 
Because of the problems with the diffractive orbits 5a', 6a' ... we ignored that series of 
orbits. Additionally, we made an approximate fix for small k of the orbits 4a, 4a' and *10b 
by ignoring 4a' and *10b and halving the amplitude of 4a. This simulates the effect of 
the cusp which has the approximate effect of halving the domain of the stationary phase 
transverse integral used to derive the geometrical trace formula. The result for the even 
spectrum is shown in Fig. ^ with and without diffractive orbits. The peaks are identified as 
corresponding to quantum eigenvalues. We observe that including diffraction helps to resolve 



the peaks but barely changes their positions. The positions are compared in Table |ITI[ The 
discrepancy is typically .06 compared to a spacing of about 0.80. Including the diffractive 
orbits changes the peak positions only by about .02 and not necessarily by the correct sign. 

This is a strange result because the amplitude of the geometric and diffractive contribu- 
tions to the density of states are comparable, as shown in the right half of Fig. ^j. We saw 
in the previous section that the eigenvalues do contain information about diffraction since 
there were diffractive peaks in the Fourier transform so it must be true that the diffractive 
orbits have some effect on the eigenvalues. To unravel this paradox it is probably best to 
calculate with zeta functions rather than traces, which we will do in a later publication |^ 



For now we note that the wedge billiard was also successfully quantised using only geometric 
orbits although it too has a diffractive vertex pB| ]. 

For completeness, we show in Fig. ^ the results for the odd spectrum and the results 
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are also enumerated in Table |nT| . Now only geometric orbits contribute. The agreement is 
somewhat better; with the exception of the third state, which is not well resolved in any 
case, the differences between the peak positions and the quantum eigenvalues are typically 
about .04. The only qualitative difference between the two parities is the diffraction so the 
difference in accuracy is presumably diffraction related. 



V. SPECTRAL STATISTICS 

It is by now well known, if not well understood, that the statistics of eigenvalues of 
chaotic systems |^2[ follow closely the predictions of random matrix theory (RMT) (for a 
review see Ref. |^3|). This has been confirmed in many examples including the conformal 
mapping of the circle of which the cardioid is a limiting case. Nevertheless, we repeat 
this here to see if the diffraction has any effect on the statistics. The fact that the diffraction 
is almost completely limited to the even spectrum means that we can study its relative effect 
by comparing the results for the even and odd spectra. A similar study was already reported 
in Ref. |^ for the half cardioid but without an explicit comparison of the odd and even 



spectra to explore possible diffraction effects. 

We first display the spacing distributions in Fig. |TT]. This is the distribution of the 
spacings s between adjacent energy levels measured in units of the local mean level density 
as found from the Weyl formula Eq. (0). The Gaussian Orthogonal Ensemble (GOE) is a 
random matrix ensemble which predicts a spacing distribution very close to 

P(s) = |sexp(-7rsV4). (38) 

Both spectra are observed to be consistent with that limit and with each other so diffraction 
appears to have no significant effect on the spacing distribution. 

The spacing distribution is a measure of short range correlations. A statistic commonly 
used to probe the long range correlations is the spectral rigidity A3(/). Here / refers to 
distances in the spectrum measured in units of the mean level spacing. (We use a small 
letter I to avoid confusion with the periodic orbit lengths which we referred to with a large 
L.) This statistic measures the average deviation of the staircase function from a local 
straight line fit over a window of length /; the average being taken as the window is moved 
through the spectrum. The GOE formula is given as a complicated integral representation 
but for large / is nearly (log/)/7r^. Chaotic systems follow the GOE result but at 



some point begin to deviate from it and finally saturate for arguments larger than about 
^max = 27r(/9(/i;))/Lj^jj^ where {p{k)) is the average level density and -Z^j^in is the length of 
the shortest periodic orbit. The shortest periodic orbit in the cardioid is 2a which has length 
2.60 in the fundamental domain. The average density of states in the range considered is 
20.0 so that /max = 48.3. 

In Fig. |T^, we show the A3(/) results for the even and odd spectra. We also show the 
saturation value of 0.265 which we found numerically. As can be seen, /max is a good 
approximation to where the saturation begins. The results are significantly different from 
the GOE prediction for values of / larger than about 7. In comparing the two data sets, 
we observe that the odd spectrum has values which are consistently larger than the even 
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spectrum for < / < 17. The difference is typically about 0.01. To determine if this is 
significant we need to compare the difference to the typical variance in A3(/). If we assume 
that each range of length / used in determine A3(/) is statistically independent, we find 
that the typical variance is 0.006. This assumption is problematic since we know that there 
are strong correlations in the spectra. If instead we determine the variance by finding the 
effect of removing selected levels, then the typical variance is 0.02 [^]. In either case, 
we conclude that if there are significant deviations between the two spectra, we are not 
sufficiently sensitive to resolve them with confidence. This conclusion is in agreement with 
the results of the wedge billiard |^ and of pseudo-integrable billiards [0]. We can increase 
the statistical significance of this statement by finding more eigenvalues. It would also be 
interesting to study this question analytically by extending the Berry's original semiclassical 
calculation of spectral rigidity to include the effect of diffractive orbits. 



VI. CONCLUSION 

In this paper, we have derived a trace formula for periodic orbits diffracted by vertices. 
The presence of the diffractive orbits causes additional structure in the quantum spectra 
of Hamiltonian systems. This structure is suppressed relative to the contributions from 
geometric orbits. The diffractive trace formula has a very similar structure to the trace 
formula for geometric orbits and from it we can find a zeta function in close analogy to 
the zeta function for geometric orbits. An important difference in the structure of these 
functions is that the diffractive zeta function involves only one product. Multiplying these 
zeta functions gives the total zeta function which will probably provide the cleanest method 
of finding the semiclassical eigenvalues. 

We specialised the discussion to the example of the cardioid billiard possessing a cusp 
at the boundary which is locally a half-plane vertex. There is overall good agreement 
between the Fourier transform of the exact spectrum and that of the trace formulae and 
we successfully reproduced geometric, diffractive and doubly diffractive peaks as well as the 
interferences among them. However, there is a region of L which is not well reproduced. Two 
reasons for this disagreement are that for certain choices of angles the diffraction picture 
breaks down and because the cusp is not a perfect half-plane vertex. 

There exists a symbolic dynamics which includes the periodic diffractive orbits in a 
natural way by inclusion of one more symbol in the alphabet. This leads to a simple result 
when discussing the symmetry reduction. For every non self-dual diffractive orbit, there 
exists a complementary self-dual one. These interfere so that the diffraction affects only the 
even spectrum and leaves the odd spectrum alone. We used the sum over periodic orbits 
to find the first few quantum eigenvalues with an accuracy of a few percent. We found the 
puzzling result that including diffractive orbits seems to have very little effect. This issue 
will be addressed in a later publication. 

In the last section we studied the level statistics of the even and odd spectra separately. 
Comparing the results from the two spectra is a probe of the effect of the diffractive orbits. 
This is because the odd spectrum has essentially no contribution from diffraction. We found 
no significant differences in either the spacing distributions or ^^{l) although we can not 
rule out the possibility that such differences might become apparent if we included more 
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states. 
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TABLES 



orbit 


ruil domain 






fundamental domain 


2a 


-1 — 






1 


4b 


+ -+- 






0101 = (01)2 


3a 


+ - + 






oil 


fib 


+ + + 






omnm — (TinTi^ 


4a 


+ + +- 






0011 


8b 


+ + + + 






00010001 = (0001)2 


8c 


+--+-++- 






10111011 = (1011)2 




TABLE I. Some geometric orbits and their 
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TABLE II. Some diffractive orbits and their 


symoois. 
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Exact 


Geometric 
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number 
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TABLE III. Top: 


even states calculated with and without diffractive orbits. 


Bottom: odd states. 
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FIGURES 



FIG. 1. (a) A path connecting the point x' to the point x via the vertex ^. The angles a, 
and Q appear in Eq. (^) defining the diffraction constant, (b) A schematic diagram of a periodic 
diffractive orbit with its local coordinates. 

FIG. 2. The cardioid billiard. A generic point ziQ") at the perimeter is shown as well as the two 
special points = and 9 = it (the cusp). We restrict 9 to the interval [— vr, vr). The dashed line 
indicates the symmetry axis. 

FIG. 3. Various geometric orbits of the cardioid billiard labeled with their names and their 
lengths. 

FIG. 4. Various diffractive orbits of the cardioid billiard labeled with their names and their 
lengths. We indicate the incoming and outgoing directions where it is ambiguous. 

FIG. 5. A family of related pruned geometric orbits and pruned diffractive orbits labeled with 
their names and their lengths. 

FIG. 6. The solid curves are the Fourier transforms of the exact spectra and the dashed curves 
are the approximations from the trace formulas Eqs. (^) and (pO|). 

FIG. 7. Comparison between the exact result (solid curve) and the diffractive trace formula 
(dashed curve) for orbit la'. The top box is the window 2.01 < k < 26.32 and the bottom box is 
the window 13.09 <k< 39.39. 

FIG. 8. As in Fig. ^ but for a range corresponding to the double diffractive orbit 3a". 

FIG. 9. Left: The solid curve is the result of using both geometric and diffractive orbits in the 
trace formula for the even spectrum while the dashed curve is from using just the geometric orbits. 
The arrows denote the exact positions of the even states. Right: A comparison of the magnitudes 
of the geometric (top) and diffractive (bottom) contributions to the density of states. 

FIG. 10. The odd density of states as constructed from just the geometric orbits and the arrows 
denote the exact positions of the odd states. 

FIG. 11. Histograms showing the spacing distributions of the even and odd spectra. The solid 
curve is the GOE result. 

FIG. 12. The A3(/) functions for the even (o) and the odd (x) spectra. The solid curve is the 
GOE result and the dashed line denotes the saturation value. 
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